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Abstract 

csf: 

The polar motion data is analyzed to obtain accurate position of the figure axis 
{S| ■ referred to the Earth-fixed frame. The variation of the figure axis should be the basic 

object to which the geophysical events are linked. By the method of rigid dynamics, 
►—5 ; the relation between the rotational and the figure axes is derived. The polar motion 

is the motion of the rotational axis on the Earth's surface and therefore the exact 

<n : 

position of the figure axis on the surface at each moment is obtained from the polar 
motion. Since the accuracy of the recent data for the polar motion is very high, the 
obtained position of the figure axis is considered to keep good quality. As an average 
for a long duration of time, the obtained figure axis exhibits stable annual and semi- 
annual variations while it has no component with the period of the Chandler wobble. 
Besides this general feature, it shows different patterns from year to year as well as 
small irregularities with shorter periods, most of which are considered to be significant 
but not attributed to the observational errors. Further, a simple model is introduced 
for the cause of the seasonal variation of the figure axis. The model explains both the 
variation of the figure axis obtained above and that of the rotational speed which is so 
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far known. 
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g ; 1 Introduction 

The polar motion is a phenomenon in the rotational motion of the Earth, appearing 
as the motion of the point on the Earth's surface at which the rotational axis pierces 
the surface. The motion is roughly circular but with changing radius. 

This phenomenon occurs in the free rotation of a rigid Earth as well. When the 
Earth is regarded as a rigid oblate spheroid, the locus is a complete circle. The period 
of the circular motion in the rigid case is A/{C — A) times the period of the Earth's 
rotation, C and A being the moments of inertia of the Earth around its symmetrical 
axis and around an axis on the equator of the symmetrical axis, respectively. We 
suppose the approximation that the Earth is a spheroid throughout the present study. 

For the value of A/(C — A) in the actual Earth, the period is about 304 sidereal 
days. The motion of the rotational axis in a rigid Earth is called the Euler motion and 
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its existence and the behavior can be proved theoretically. 

In the actual Earth, the Euler motion suffers a considerable deformation. The 
motion is now not with a simple shape and with a single period but a compound of 
the various elliptical motions with different periods, shapes and sizes, all varying with 
time. 

The largest ellipse among them stems from the Euler motion in the rigid case and 
is called the free oscillation (Munk and McDonald, 1960). The period of of the free 
oscillation is about 435 days, considerably longer than that in the rigid Earth. The 
free oscillation is also called the Chandler wobble after the name of the discoverer of 
the period. 

The prolongation of the period in the free oscillation in the actual Earth is attributed 
to the elasticity of the Earth. It was first explained by Newcomb (1892). Explanations 
are also seen in Munk and McDonald (1960) as well as in Kubo (1991). As it is 
important to take the prolongation of the period into account in the argument in 
Section 2, a detailed review of it is given there. 

Besides the Chandler wobble and the many other periodical components in the 
polar motion is also found a secular drift of its center. 

All the motions other than the free oscillation are considered to be reflecting the 
variation of the mass distribution in the whole Earth including the ocean and the 
atmosphere. Therefore, the investigation of the polar motion has held a very important 
position for understanding the dynamical state of the Earth. 

For that reason, numerous analyses have been published concerning the polar mo- 
tion, especially for its periodic parts. Among them, some studies aim only to analyze 
the pattern of the polar motion itself such as the variations in its amplitude and phase, 
e.g., Kimura (1917), Guinot (1972), Hopfner (2002). Some other studies intend to ex- 
plain the variation in the polar motion by geophysical causes such as the atmosphere 
(Aoyama and Naito, 2001), the ocean (Dickman, 1993; Gross, 2000), geomagnetic jerks 
(Gibert, 2008) and so on. 

The polar motion, however, is considered to reflect the state of the Earth not so 
directly as the variation of the figure axis. The change in the inner state of the Earth 
first brings about a variation in the figure axis, which then affects the polar motion. 
Therefore it is the variation of the figure axis on the Earth's surface but not the polar 
motion that geophysical events should be linked to. 

The present study intends to obtain the position and the motion of the figure axis 
relative to the Earth's surface. It is performed by the method of rigid dynamics, 
regarding the elasticity and anelasticity of the Earth as perturbations to the rotation 
of a rigid Earth. 

Primarily, a main theme of the rotation of the Earth in rigid dynamics is to obtain 
the respective motions of the rotational axis and the figure axis. The solutions to 
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the rotation of the rigid Earth as well as of the elastic Earth have been published 
numerously so far, among which a typical example for the former is Kinoshita (1977) 
and one for the latter would be Kubo (1991). 

It is not necessary at present, however, to solve the respective motions of both axes. 
We only have to determine the position of the figure axis at any moment on the Earth's 
surface and we can attain this purpose only if the relation between the rotational and 
the figure axes is known because the observed polar motion is the position of the 
rotational axis on the Earth's surface. Both the figure axis and the rotational axis are 
definite physical entities and the relation between them can be determined without any 
ambiguity. 

The derivation of the fundamental equations to give the relation between both axes 
is carried out in Section 2 based on the result in Kubo (1991), and some preliminary 
studies using the equations are presented in Sections 3 and 4. 

In Sections 5 and 6, the actual polar motion of the Earth is exermined. First in 
Section 5 are overviewed the feature of the polar motion published by IERS as well as 
the spectrum distribution for the components contained in it, with a result that all the 
conclusions reached are not different largely from those preceding. 

In Section 6, on the other hand, is presented a quite new result on the position and 
the motion of the figure axis which are obtained from the data of the polar motion. 

Finally, in Section 7, a simple model of the mass transfer in the Earth is introduced 
to examine if the variation of the figure axis obtained above is reasonable or not. 

2 Equations for the polar motion 

2.1 Definitions of the figure axis and the rotational axis in the 
present study 

The polar motion in the rotation of the Earth is defined as the motion of the rotational 
axis referred to a coordinate system set on the Earth surface near the north pole. We 
separate this motion into two parts: One is the motion of the Earth's figure axis in 
the said coordinate system and the other is the motion of the rotational axis relative 
to the moving figure axis. The motion of the figure axis does not exist in the rotation 
of the rigid Earth but only the latter motion. 

In the deformable Earth, where the variation of the figure axis appears, it is impor- 
tant to give exact definitions for both the figure axis and the rotational axis. Therefore, 
first we consider the rigorous definition for both axes adopted in the present study. 

As for the rotational axis in the present study we adopt Celestial Intermediate Pole 
(CIP) which is defined definitely in the resolution of the International Astronomical 
Union (IAU, 2000). 
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CIP, which was called Celestial Ephemeris Pole (CEP) before, defines the fundamen- 
tal coordinate system of date in the celestial sphere. The precession and the nutation 
in the astronomical ephemeris are defined as the secular and the periodic motions of 
CIP, respectively. CIP is the axis around which the Earth actually makes the diurnal 
motion but different from the rotational axis in mathematical sense by less than 0.01" 
(Seidelmann, 1982; IERS, 2004). 

The motion of the rotational axis defined by CIP does not contain diurnal compo- 
nent in space nor with respect to the surface of the Earth. Therefore the polar motion 
defined by this rotational axis does not have diurnal component. The coordinates of 
the pole in the polar motion published from IERS are based on this axis. 

Next we define the figure axis in the present study. The figure axis of the Earth 
which is changing at every instant is the axis with the largest moment of inertia at this 
instant. In the present study, however, we do not consider changes with short period 
close to the diurnal motion, for example due to the body tide and others (if they exist) 
but the axis averaged as for the periodic changes, say, shorter than several days. 

More importantly, however, we must take into account the fact that the actual 
figure axis at any instant is the axis with the largest moment of inertia formed as the 
compound of the Earth's figure which the Earth would have if it were not rotating, i.e., 
the intrinsic figure of the Earth, and the bulge produced by the centrifugal force due to 
the rotation which is formed symmetrically around the rotational axis (the rotational 
axis defined just above). In the following we call this axis for the compounded figure 
of the Earth the effective figure axis. 

The effective axis displaces from the intrinsic axis by an amount explained later. 
This displacement does not occur with a short period but a comparably long period, 
synchronous with the period of the Chandler wobble. However, we do not let this 
displacement of the effective figure axis be included in the definition of the figure axis 
in the present study. That is, we define the figure axis by the intrinsic figure axis 
but not by the effective figure axis. This definition of the figure axis is essentially not 
different from the generally accepted one. 

Anyway, if we know the location of the figure axis thus defined, that of the effective 
figure axis can be found easily by following the procedure inferred by the description 
in the next subsection. 

2.2 Equations describing the rotational axis relative to the 
figure axis 

Now we consider the relation between the rotational axis and the figure axis defined in 
the previous subsection. Before discussing the subject, however, it will be significant 
to glance over the nature of the polar motion in the rigid Earth as well as in the 
deformable Earth which the rigid dynamics describes. 
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Figure 1: The Earth's surface near the north pole. C and E are the intrinsic and 
the effective figure axes at some instant, respectively. P is the rotational axes at the 
same moment and P' is that some time duration before. The origin of xy-coordinate 
system is at IRP (International Reference Pole) and x- and y-axes in the directions of 
Greenwich and 90° east longitude, respectively. 

As mentioned in Section 1, in the free rotation of the rigid Earth, the point on 

the Earth's surface at which the rotational axis pierces the surface rotates around the 

C — A 

figure axis with the angular speed — — — a;, uj being the sidereal angular speed of the 
Earth's rotation (Kinoshita, 1977). This circular motion of the rotational axis on the 
surface of the rigid Earth is nothing else but the Euler motion. 

A similar circular motion of the rotational axis occurs also in the deformable Earth. 
Suppose that the Earth is deformed by the centrifugal force due to the rotation and 
then consider a rigid Earth with the same values of A and C as in the deformed Earth 
and its figure axis coincident with the effective figure axis, i.e., a rigid Earth equivalent 
to the deformed Earth. Then the rotational axis of the deformed Earth behaves in the 
same way as in the equivalent rigid Earth at each instant (Kubo, 2009). 

To return to the subject, we consider first the Earth which has elasticity but not 
plasticity. Fig. [I] shows the Earth's surface near the north pole. Let C and P be the 
figure axis and the rotational axis at some instant and let CP = j. The bulge by the 
centrifugal force due to the rotation is produced at the equatorial zone of the rotational 
axis. Acoording to Kubo (1991), the figure axis of the Earth produced as the conpound 
of the intrinsic Earth's figure and the bulge due to the centifugal force, i.e., the effective 
figure axis, is located at the point E on the line CP and at the distance ej from C with 

3/iC 

where \i is given by 

" = = °- ooii6fc ' (2) 
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k being the Love number and a e the equatorial radius of the Earth, respectively, and 
G the universal gravitational constant. 

Then, as mentioned above concerning the behavior of the rotational axis toward 
the effective figure axis in the deformable Earth, P moves so as to rotate around E 
with the angular speed 

O — A 

Q = ^—r^OJ. (3) 

A v J 

It should be noticed that A and C in this equation is those for the actual Earth shape 
including the centrifugal effect, that is, the observed values of the moments of inertia. 
Therefore the value of Q is (l/304)u;, or the corresponding rotational period is about 
304 sidereal days. This is the rotational speed of the Euler motion in the rotation of 
the equivalent rigid Earth. 

As a consequence that P makes sush a motion, P rotates around C with an angular 
speed 

Q! — (1 — e)Q = (1 ~ 3 ^ C ~^ . (4) 

Supposing k = 0.29, this angular speed corresponds to the period of about 440 days, 
the period of the Chandler wobble. Thus the period of the Euler motion is prolonged 
to that of the Chandler wobble due to the elasticity of the Earth, being coincident with 
the expression in Kubo (1991) which was derived by a more analytical method. 

Next we consider the Earth which has anelasticity besides elasticity. We treat this 
condition by introducing a model such that the bulge by the centrifugal force at a 
moment is produced around P', the rotational axis some time interval ago, but not 
around P, the rotational axis at the moment. Let ZPCP' = S. Then the effective figure 
axis at the moment comes to the point E' on the line CP 7 and at a distance ej from C. 
We consider that and therefore ej as well, does not change during the time interval 
S/uj corresponding to the motion from P' to P, which is probably one day or so. This 
results in that P moves so as to draw a circle around E' with the radius E'P, which is 
equal to (1 — e)j with a sufficient approximation, and with the angular speed Q. 

This situation is expressed by the following equations: We consider a coordinate 
system on the Earth's surface in Fig. [TJ Although this coordinate system is arbitrary, it 
may be convenient to choose International Reference Pole (IRP), which can be regarded 
as practically the same as Conventional International Origin (CIO), as the origin and 
to take x-axis in the longitude of Greenwich. Then y-axis is in 90° east longitude. 

Let the coordinates of P, C and E' be (x, y) ) (x c , y c ) and (x e , y e ), respectively. Then 
the velocity of P relative to E' is written as 

i = -{y-y e )n, /rX 
y = (x- x e )Q. 

We then express (x e) y e ) in terms of (x C) y c ) in these equations. Let the angle of the 
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direction of the line CP from x-axis be 9. Then 



x e = x c + ej cos(# — 5) 




(6) 



= x c + e[(x - x c ) + (y- y c )S]. 



Therefore we have with a sufficient accuracy 



x - x e = x - x c - e[(x - x c ) + (y- y c )S] 
= (1 - e)(x - x c ) - e(y - y c )S. 



(7) 



Similarly, 



y e = (1 - e)(y - y c ) + e(x - x c )5. 



(8) 



Hence, 



x 



y 



-[(1 -e)(y- y c ) + e(x - x c )5]Q, 
[(1 - e)(x - x c ) - e(y - y c )5]Cl. 



(9) 



These are the equations for the polar motion which gives the motion of the rotational 
axis P(x, y) when the location of the figure axis C(x ci y c ) is given as a known function 
of time in the Earth with the elasticity and anelasticity represented by /i and 5, re- 
spectively, or inversely gives the position of the figure axis when the polar motion is 
known. 



Now, as a preliminary study to use the equations derived above, we try to generate 
a polar motion by means of numerical integration giving an arbitrary motion of the 
figure axis. 

For example, we give a motion of the figure axis as follows: 



where v is the angular speed corresponding to the period of one tropical year, i.e., 
27r/365.2422 days. f(t) is a step function which takes the values of either +1 or — 1 
changing them with random durations of time of 50 x i days (i = 1 to 9). The motion 
of the figure axis f fTOl ) is shown in Fig. [2l 

The variations in x c is to give an annual variation in the generated polar motion 
as seen in the actual polar motion while f(t) in y c is to give an irregularity to the 
generated polar motion. Without this function, the generated polar motion would be 



3 Artificial polar motion 



x c = O.OS^cosz/i, 
y c = 0.01" x f(t) 



(10) 
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Figure 2: The coordinates of the figure axis arbitrarily given for 20 years, x-coordinate 
represents the annually periodic nature of the variation of the axis and y-coordinate 
an irregularity in it. 

a simple periodic motion when 5 = and would converge to a periodic motion even if 
5 Cu. The annual term and the irregular term in the motion of the figure axis are 
separated into x c and y c in order to make it easy to conferm how faithfully the figure 
axis is reproduced from the polar motion in the next section. 

As for the initial coordinates of the rotational axis P, we set them as 

x = +0.1", y = 0.0". (11) 

These initial values are arbitrary but not given so that the generated polar motion fits 
the actual one in the Earth, that is, the generated polar motion is quite fictitious. 

We also give the Love number 0.29 which corresponds to the values of \i — 0.000336 
and e = 0.308 in performing numerical integration. The corresponding period of Chan- 
dler wobble is about 440days. As for the value of S which is a parameter related to the 
damping coefficient of the Earth, we first adopt it as 0. 

A numerical integration by the classical Runge-Kutta method for the time span of 
20 years (7305 days) gives a polar motion, i.e., a variation of the coordinates of the 
rotational axis (x, y) as shown in Fig. [3j The power spectra for the variation are also 
shown in the figure. 

We next carry out the same numerical integration but changing the value of 5 to 
1.0°. The variation of (x, y) and their spectra in this case are shown in Fig. 0J 

From the result of the numerical integration above we see a difference between the 
polar motions obtained assuming 8 = and S = 1.0°, but it is difficult to judge which 

1 If x c and y c take the form of x c = a smut + bcosvt and y c — csmvt + dcosvt in Eq. (j9]), the 
equations have the following general solution: 

x = de~ emt cos[(l - e)Qt + C 2 ] +psmvt + q cos ist, 
y = de~ emt sin[(l - e)Qt + C 2 ] + r sin ut + s cos vt, 

where C\ and C 2 are the constants of integration and p, r and s are linear combinations of a, 6, c 
and d with the coefficients being functions of e, 5, Q and v. If 5 = 0, obviously x and y are periodic 
functions. If S ^ (and S > 0), the exponential terms in the solution attenuate and x and y converge 
to periodic functions, too. 
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Figure 3: The time variation of the coordinates and its power spectra for an artificially 
produced polar motion with the figure axis given in Fig. [2] and assuming 5 = 0. 





X 












y 









0.3"- 




0.2 - 




0.1 - 









-o.i - 


CD 


-0.2 - 








-0.3 






< 


0.3" 




0.2 - 




0.1 - 




z 




-o.i - 




-0.2 - 




-0.3 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 

Time 



19 year 




100 150 



250 300 350 400 450 500 550 
Period 



Figure 4: The time variation of the coordinates and its power spectra for an artificially 
produced polar motion with the figure axis given in Fig. [2] and assuming S = 1°. 
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represents the actual polar motion better from this result. Due to the irregular function 
in ( fTUl) . the free oscillation does not attenuate in both cases. In other words, it would 
be difficult to determine the value of 5 from the data of the actual polar motion. 

We further notice in the spectra both in Figs. [3] and [4] that the period of the free 
oscillation appearing in the generated polar motion is close to 440 days, which is the 
same as the assumed value in the numerical integration. This fact is not necessarily 
self-evident. Although it may be possible to prove it theoretically anyhow, we can say 
that we have shown it empirically here. 

4 To reproduce the figure axis from the generated 
polar motion 

Next we consider an inverse problem of the one in the previous section. That is, 
contrary to the process to generate a polar motion giving a motion of the figure axis, 
we consider the problem to obtain the location of the figure axis at each instant from 
the produced polar motion. 

Fig. [5] shows the same area as Fig. [U Let C(x c , y c ) represent the figure axis at an 
instant t\ and Pi(xi,?/i) and V^ix^V'i) the rotational axes at t\ and t 2 ? h being the 
instant t days later than t\. We approximate that the figure axis does not move during 
the period of t days and the rotational axis draws a quasi circle around the figure axis 
decreasing the radius gradually due to the attenuatior§. 

Define angles £,77 and £ as in Fig El Assume T (days) which is the period of the 
circular motion of the rotational axis P around the figure axis C corresponding to the 
angular speed of P at the instant, and let P1P2 be d. Then ( is approximately expressed 
by 

2nd , , ... 
C = — (radian). (12) 

Then, from the cosine formula for the plane triangle CP1P2, 

d 2 =f 1 +f 1 e-^-2f 1 e-^cosC, (13) 



ji being CPi. From this 



Ji = = • (14) 

^l + e -2at _2e- ff *cosC 



Applying the sine formula to the same triangle: 

sin i] sin ( 

jie- at ~ d ' 



(15) 



2 If CPi = j, then CP2 = je at , a being the damping coefficient denned by —j/j. According to 
Kubo (1991), a = 1.12 x 10 -4 £/day, S being expressed in degree of angle. S = 1° corresponds to the 



damping with the half-life of about 24.5 years. For t = 10 days, for example, CP2 = 0.9989 x CPi. 
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Figure 5: The same area as in Fig. [TJ Pi and P 2 are the rotational axes at time t\ and 
h (h < t 2 ). The figure axis C is supposed to be constant during the time duration t\ 
to t 2 . 

we have rj. Meanwhile £ is obtaied by 

tan^ = (16) 
X\ — x 2 

Usind r] and £ thus obtained, we have the coodinates of C by 

%c = Xl ~ jl C0S(£ + T?), 

Vc = Vi ~ji sin(f + 77). 
These coordinates are regarded as those of the figure axis at the time — ~^ 



(17) 



2 

Now, following this procedure, we obtain the locations of the figure axis for every 
day from the data of the artificial poler motion generated in the previous section. In 
doing this, of the two polar motions we obtained with different values of 5, we adopt 
the one given in Fig. 0J which is generated with the value of S = 1.0°. 

Further, in order to perform this, it is necessary to assume the period of the free 
oscillation T in Eq. (|T2l) . Redgarding this, we may assume 440 days which has been 
obtained in the spectrum analysis (Fig. H]). As seen in the previous section, the value 
appearing in the spectrum shows a very good coincidence with the value given for the 
period of the free oscillation when producing the artificial polar motion. 

As for the damping coefficient 8 (therefore for a), on the other hand, we can not 
know the true value for the actual Earth. So, concerning the artificial polar motion, 
we analyze it first assuming the value of and then assuming the value 1°, the latter 
being the same as the one used in the simulation, i.e., the true value so to speak. 

The result is shown in Figs. [6] and [3 for 8 = and 5=1°, respectively. In the two 
figures, Figs. [6]and[7l we do not see a large difference between the cases assuming 5 = 
and 8 = 1.0°, and the original motion of the figure axis given by Eq. ([2]) is reproduced 
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Figure 6: The coordinates of the reproduced fgure axis from the polar motion shown 
in Fig. [Hand assuming 5 = 0. 
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Figure 7: The coordinates of the reproduced fgure axis from the polar motion shown 
in Fig. [Hand assuming 5=1°. 

very well in both cases, although we admit a slight superiority of the reproduction in 
the case where the true value is assumed for the damping parameter. This suggests 
that the choice for the value of the damping parameter which must be assumed in the 
process of obtaining the figure axis from the data of the polar motion does not affect 
the final result so much. 

Anyway, it is shown through this example that the figure axis can be determined 
from the polar motion with a very good accuracy. 

5 The polar motion in the actual Earth and its 
primitive analysis 

In the following sections we treat the actual polar motion observed in the Earth. The 
data of the polar motion or the coordinates of the rotational axis is acquired as Bulletin 
B data published by the International Earth Rotation Service (IERS), which can be 
downloaded from the internet website (IERS, 2012). Of these data we use the daily 
data for the period of 20 years (7305 days) from 1 January 1992 to 31 December 2011. 
It is emphasized that they are all compiled based on the modern observation data 
acquired using the space-geodetic techniques. 

First we look at the original data of the coordinates of the rotational axis (x, y), 
which is shown in Fig. El The origin of the coordinate system is the International 
Reference Pole (IRP) and x- and y-axes are taken in the directions of the Greenwich 
longitude and 90° east longitude. So, attention should be paid to that the sign of 
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Figure 8: The time variation of the coordinates and its power spectra for the actual 
polar motion observed in the Earth for the recent 20 years. 

y-coordinate is opposite to that in the IERS original data. The power spectra of x- 
and y-coordinates are shown also in Fig. [8] 

The polar motion consists of a secular motion and periodic motions. Concerning 
the periodic variation, they consist of mostly the free oscillation (the Chandler wobble), 
the annual and the semi-annual components, among which the semi-annual component 
is comparatively small. 

Looking at the spectra, we notice first the prevailing Chandler period of about 437 
days. And, other than this main component, the annual and semi-annual components 
are recognized. On the whole all these features are consistent with many preceding 
analyses, although some minor differences may exist due to the difference of the period 
for the analysis and so on. 

Regarding the coordinates (x, y) of the IERS data shown in Fig. [8j we express 
them as functions of time assuming the quadratic polynomial for the secular motion 
and fitting respective ellipses for the periodic motions of the Chandler wobble and the 
annual and semi-annual terms. That is, we perform the least squares fitting assuming 



13 




Figure 9: The Chandler (left) and the annual (right) components in the polar motion 
averaged for 20 years, the semi-annual component being incorporated in the latter. The 
0-th day in the Chandler component corresponds to 1 January 1992 and the days after 
every Chandler period. The names of the months in the annual component represent 
the first day of each month. 



the expressions for x and y, respectively, as 

x = a + ait + a 2 t 2 
+ C\ sin Q!t + c 2 cos Q!t + c 3 sin vt + C4 cos vt + c 5 sin 2vt + c 6 cos 2z/£, 

(18) 

y = b + M + 62* 
+ c?i sin fl't + d,2 cos Q't + o?3 sin z/t + <i 4 cos i^t + <is sin 2i/t + <i 6 cos 2vt, 

where t is the time beginning from 1 January 1992, expressed in one tropical year, i.e., 
365.2422 days. Also Of = (2vr x 365.2422/437) radian/year and v = 2ir radian/year. 

The result is: 
For the secular motion, 

a = [+0.034(±0.002)f, b = [-0.318(±0.002)]", 

ai = [+0.0022(±0.0008)]"/t, 6i = [-0.0023(±0.0008)] ,/ /t, (19) 

a 2 = [-0.00036(±0.00009)] ,/ /t 2 , b 2 = [-0.00013(±0.00009)] ,/ /t 2 . 

For the Chandler wobble with the period of 437 days, 

ci = [-0.1001(±0.0007)]", di = [+0.1036(±0.0007)]", 
c 2 = [+0.1039(±0.0007)f, d 2 = [+0.0986(±0.0007)]". 

For the annual variation, 

c 3 = [-0.0777(±0.0007)f, d 3 = [-0.0424(±0.0007)]", 
c 4 = [-0.0462(±0.0007)] // , d 4 = [+0.0722(±0.0007)f. 

For the semi-annual variation, 

c 5 = [+0.0043(±0.0007)]", d b = [+0.0005(±0.0007)]", 
c 6 = [-0.0009(±0.0007)] // , d 6 = [-0.0014(±0.0007)]". 



(20) 



(21) 



(22) 
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Eqs. f |20l) to f |22l ) give the locii of the free oscillation and the annual and the semi- 
annual terms and they are shown in Fig. [9j in which the semi-annual component is 
included in the annual one. The two locii exhibit nearly circular patterns. 

Here, notice should be taken that the patterns are the averaged ones for the duration 
of 20 years from 1992 to 2011. That is, the periodic components are assumed roughly 
constant throughout the period. As for the secular motion we consider the problem 
below again. 

6 The position and the motion of the figure axis 
obtained from the polar motion 

Now we apply the procedure described in Section 4 to the actual polar motion shown 
in Fig. [8] and obtain the positions of the figure axis on the Earth's surface for every 
day during 20 years. As for the period of the free oscillation in the polar motion we 
assume 437 days, the same as obtained in the spectra in Fig. [8J since we know that the 
period of the free oscillation to assume in Eq. ( JT2l) appears as the peak in the spectra 
as we saw in Section 4. 

Also we do not know the value for S. As for it we assume 0. Concerning this we 
saw also in Section 4 that the error in this choice does not affect the result so much. 

The obtained values of the coordinates x c and y c for every day and the spectra 
for the periods contained in their variations are shown in Fig. [TOl Remarkably, there 
exists no component with the period of the free oscillation around 437 days at all in 
the spectra. This should be natural because it is not likely that the variation of the 
figure axis is affected by the free oscillation. Therefore, in finding the expressions as 
functions of time for those coordinates by carrying out the least squares fitting, we 
assume forms similar to Eq. f fT8l) but not containing the oscillation with the period of 
437 days, i.e., not containing Ci,c 2 ,rfi and d 2 \ 

%c — a + a lt + a 2t 2 

+ c 3 sin vt + C4 cos vt + c 5 sin 2vt + c 6 cos 2vt, 

(23) 

y c = b + b x t + b 2 t 

+ d 3 sin vt + rf 4 cos vt + d 5 sin 2vt + d 6 cos 2vt, 

The result is: 
For the secular motion, 

a = [+0.035(±0.001)] // , b = [-0.327(±0.001)] // , 

ai = [+0.0004(±0.0003)]"/t, b x = [+0.0011(±0.0006)] // /t, (24) 

a 2 = [-0.00001(±0.00004)] // /t 2 , b 2 = [-0.00045(±0.00007)] // /t 2 . 
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Figure 10: The time variation of the coordinates and its power spectra for the figure 
axis obtained from the polar motion shown in Fig. [8] and assuming 5 = 0. 

For the annual variation, 

c 3 = [+0.0085(±0.0003)] 
c 4 = [+0.0004(±0.0003)] 

For the semi-annual variation, 

c 5 = [-0.0001(±0.0003)] 
c 6 = [-0.0046(±0.0003)] 

In examining the result in the analysis in this section in reference to that in the 
previous section, we first recognize an innegligible difference of about 0.01" between 
the constant terms in (x,y) and (x c ,y c ) given by Eqs. (|T9l) and ( 12^1) . respectively. 
Concerning that, we must say that the averaged position of (x C) y c ) has a significant 
physical meaning but that of (x, y) does not necessarily. As for the secular motions, 
on the other hand, it will be difficult to say something significant. 

Next, when we compare the annual and the semi-annual variations in (x,y) and 
{x C) y c ) ) we see a striking difference. The amplitudes of the annual and semi-annual 
components in (x C) y c ) given by Eqs. ( |25l) and ( |26l) are much smaller than those in (x, y) 
given by Eqs. ( |2B and ( |22l) . Further, the locus of the annual and the semi-annual 
components given by Eqs. ( |25l) and (|26l) and shown in Fig. [TT] is not like a circle but 
rather close to a line compared to that in the polar motion. 



d 3 = [+0.0102(±0.0005)]", 
d 4 = [-0.0236(±0.0005)] // . 



d 5 = [+0.0038(±0.0005)] // , 

d 6 = [+0.0065(±0.0005)] // . >j 
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Figure 11: The locus for the motion of the figure axis in the Earth, averaged for 20 
years. The names of the months represent the first day of each month. 

Here again it should be noticed that the locus is the averaged one for 20 year 
duration. The shape changes from year to year and, other than that, there exist many 
small variations with shorter periods. 

Since the modern observation of the polar motion has an accuracy as high as 0.0001" 
and the time resolution as short as one day, the position of the figure axis obtained from 
the polar motion data is considered to have a corresponding good accuracy. Although 
only the variation of the figure axis in a global time scale has been investigated here, 
the most of the details in the variation should not be attributed to the errors in the 
observations or in the process of the analysis but should be considered to reflect the 
dynamical state of the Earth. Therefore they are necessary to be examined more fully 
in future studies. 



7 A model to explain the variation of the figure 
axis 

In this section we propose a model of the Earth which brings about the variation of 
the figure axis obtained in the previous section and shown in Fig. [TU In doing this, 
it would be a touchstone to judge the validity of the model if it produces a reasonable 
magnitude as well as a rather linear locus regarding the variation of the axis. 

We now consider a case where a seasonal transfer of matter in the vertical direction 
occurs in some region on the Earth's surface. Most of the matter may be water and 
the transfer may take place either in plants or in the air, or in both. 

Suppose that the total mass to move is m and the transfer occurs in some area 
around a point with the longitude and latitude (A, 0), and let the vertical distance of 
the transfer h(t) from the mean height be expressed by 

h(t) = /sin(z/t + (27) 

if) being an arbitrary phase angle. 
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With respect to a geocentric rectangular coordinate system with x-, y-axes in the di- 
rections of the Greenwich longitude and 90° east longitude on the equator, respectively, 
the coordinates of the point are 

x = a cos (j) cos A, y — a cos sin A, z = a sin 0, (28) 

and the components of the transfer vector of the matter, Ax, Ay and Az ) are given by 

Ax = h(t) cos (j) cos A, Ay = h(t) cos (j) sin A, Az = h(t) sin 0, (29) 

where a is the Earth's radius at the point but it may be regarded as equal to the 
equatorial radius a e . 

If the moments of inertia for x-, y- and z-axes be respectively A, B and C, then the 
changes of them due to this transfer of the matter are 

AA = Am(y 2 + z 2 ) = 2m(yAy + zAz) = 2mah(t)(l - cos 2 0cos 2 A), 

AB = Am(z 2 + x 2 ) = 2m(zAz + xAx) = 2mah(t)(l - cos 2 0sin 2 A), (30) 

AC Am(x 2 + y 2 ) 2m(xAx + yAy) 2mah(t)cos 2 (f). 

At the same time, the products of inertia for x-, y- and z-axes, respectively denoted by 
D, E and F, which are originally 0, are generated as follows; 

D —A(myz) —m(yAz + zAy) —2mah(t) sin (j) cos (j) sin A, 

E —A{mzx) —m{zAx + xAz) —2mah(t) sin cos cos A, (31) 

F — —A(mxy) —m(xAy + yAx) — 2mah(t)cos 2 (f)sm A cos A. 

We next investigate how the location of the figure axis in the Earth shifts accompa- 
nying to this change of the moments and the products of inertia. We notice that this 
change of the location of the figure axis is similar to that caused by the deformation 
in the elastic Earth due to the body tide originating from the attraction of the outer 
bodies, studied by Kubo (2009) in detail. So we follow this preceding study in the 
following. 

In the non-deformed Earth or the avarage shape of ther Earth, let the axes which are 
associated with the moments of inertia A, B and C be denoted by A-, B- and C-axes, 
respectively, and those of the deformed Earth by A 7 -, B'- and Coaxes, respectively. 

The shifted new figure axis, i.e., C'-axis, is expected to be located on the meridian 
plane of the longitude A, and one of the other new principal axes, B'-axis, is also on this 
plane and separated from C'-axis by 90°. Then the other new principal axis, A 7 - axis, 
is on the original equator and in the longitude A — 90° (Fig. [T2l) . 

Consider a new coordinate system with x 7 -, y'- and z'-axes coinciding with A 7 -, B 7 - 
and Coaxes, respectively. 

Let the co-latitude of z'-axis in xyz-system be 7. The values of the moments and 
the products of inertia with respect to x'yV coordinate system, A \ £>', C", D', E f and 
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Figure 12: The averaged (ABC) and the deformed (A'B'C) figures of the Earth. The 
figure is represented by A-, B- and C-axes which are associated with the moments of 
inertia, A, B and C, respectively. 7 is the shift of C-axis due to the deformation. It 
is so assumed that the vertical transfer of natter around the point (A, (j)) causes the 
deformation. 

F 1 are obtained by examining how the equation for the quadratic surface (the inertia 
ellipsoid) : 



(A + AA)x 2 + (A + AB)y 2 + (C + AC)z 2 + 2Dyz + 2Ezx + 2Fxy = 1 
is changed to the new equation of the quadratic surface: 



A'x z + B'y u + C'z' z + 2D'y'z f + 2E , z , x f + 2F'x'y' = 1, 
after the rotation of the coordinates sytem from xyz to x'yV given by 

/ x \ ( x' \ 



y 



= Rz (g ° - A )i4( 7 ) 



V 



where Ri(£) means the rotation around i-axis by the angle £. 
Explicitly SMty is written as 



/ x \ /sin Ax' + cos A cos jy' + cos A sin 72' ^ 

y = — cos \x' + sin A cos jy' + sin A sin jz' 
\ z J y — sin jy' + cos jz' 

( sin Ax' + cos Xy' + 7 cos Xz' \ 
= — cos Ax' + sin Xy' + 7 sin Xz' 

\ -IV' + z' J 



J 

+ 0( 7 2 ). 



(32) 



(33) 



(34) 



(35) 
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Then, taking into account that 7 is the same order as mah(t)/C, we have after 
some manipulation 

A' = A + 2mah(t) + 0(7 2 ), (36) 

B' = A + 2mah(t)sm 2 (p + 0(j 2 ), (37) 

C = C + 2mah(t)cos 2 <j) + 0( 7 2 ), (38) 

D' = 2^(A -C)- 4mah(t) sin cos + 0(t 2 ), (39) 

E' = 0( 7 2 ), (40) 

F' = 0( 7 2 ). (41) 

In (1391) . if 7 satisfies the following relation: 

2mah(t) . . . 
7 = sin cos 0, (42) 

then D' = 0(7 2 ) and therefore x'-, y'- and z'-axes are the new principal axes of inertia 
to the accuracy of 0(7). 

According to the discussion above, we can conclude that the position of the figure 
axis shifts with respect to its averaged position by the angular distance 7 and in the 
direction of the longitude A. 

TYlil 

We suppose, e.g., that = 2 x 10 -10 , where m e is the mass of the Earth, and 

(j) = —15°. Then, taking into account C = 0.33m e ag, we have 

7 ^ 0.02"sin(z4 + 7/0, (43) 

almost the same magnitude of the variation in the figure axis we obtained in the 
previous section. 

If we assume I = 6m, supposing the case of the mass transfer in plants, then 

777/ 777/ 

— = 2 x 10 -4 and if I = 6000m, supposing the transfer in the air, then — = 2 x 10 -7 . 

Of course, the transfer of the mass does not occur only at one place but at many 
places. In this sense m may be written as m symbolically. And actually, it would be 
occurring in all the places on the Earth's surface and there would be yielded an averaged 
transfer which does not depend on the longitude, too. Then each local transfer which 
brings about the variation of the figure axis could be regarded as a deviation from the 
averaged transfer. 

The actual variation of the figure axis is the vector sum of many variations each 
caused by a local deviation from the averaged transfer. The shape of such a sum would 
be a line if all the local transfers are synchronous, i.e., if ip in Eq. f |27l ) is equal at all 
points, but generally be an ellipse with considerable eccentricity. 

On the other hand, a global vertical transfer of the mass independent of longitude 
causes a change of C and consequently a change of the Earth's ratational speed. For 
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simplicity, consider a case where such a global transfer occurs uniformly in the zone 
along the equator with the mass per unit area m! and the total mass M, and the 
amount of the transfer follows the same form as Eq. f l27l ). Then, very roughly, AC(t), 
the variation of C due to this transfer of mass, is given by 

AC = J 2rriah(t)d</>d\ = 2ah(t) J m'd(t)d\ = 2Mah(t). (44) 

The consequent variation of the rotational speed Auj is 

» H H Ari AC 2Mah{t) 2MaV . , in 

= ~c = ~c> = ~~c~ u = c~ 00 = — c~ u sm( + ^' ( ^ 

where H is the total angular momentum of the rotation which is constant and uj is the 
rotational speed. 

On the other hand, the seasonal variation of the rotational angle, or that of the 
universal time, is known to be comparatively stable and the following expression has 
been used for a considerably long time (USNO and HMNAO, 2000): 

AT S = 22 sin vt — 12 cos vt — 6 sin 2vt + 7 cos 2vt (ms) . (46) 

Differentiating this equation we have for the variation of the rotational speed only 
with regard to the annual terms, 

Auj = (4.4 x 10" 9 cos vt + 2.4 x 10" 9 sin ist)uj + . . . 

(47) 

= 5.0 x 10- 9 cjsin(z4 + 61°) + . . . . 

Comparing this with Eq. f l45l ). we see that they coincide roughly when we assume 

= 8 x 10 -10 and ?// = 241°, where we have taken into acount C = 0.33m e ao again. 

m e a 

Thus, assuming V = /, the ratio of the masses involved in the local deviations of the 
transfer to that in the global transfer is 

^— ^ 0.25. (48) 
M v 7 

This value of the ratio is reasonable. Therefore it is natural to conclude that the 
variation of AT S is caused by the globally averaged transfer of mass while the motion 
of the figure axis by the local deviations of the transfer. 



8 Conclusions 

The variation of the figure axis reflects the physical state in the Earth more directly 
than the polar motion since the change inside the Earth first brings about a variation 
of the figure axis and then it affects the polar motion. 

In the present study, the equatons to give the relation between the figure axis and 
the rotational axis in the Earth with both elasticity and anelasticity have been derived. 
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Since the polar motion is the motion of the rotational axis on the Earth's surface, the 
obtained relation makes it possible to determine the position of the figure axis in the 
Earth-fixed frame at any moment from the data of the polar motion. 

The position and the motion of the figure axis obtained from the actual data of the 
polar motion exhibits various characteristics. 

First, the averaged position of the polar motion are not coincident with that of the 
figure axis exactly but is dislocated by about 0.01", while the differnce between the 
secular motions of both is not so definite. 

The periodic variation of the figure axis has obvious and comparatively stable com- 
ponents of annual and semi-annual periods but no component with the period of the 
free oscillation. 

On the whole, the annual and the semi-annual components draw a shape resembling 
a thin ellipse with the semi-major axis of about 0.02". 

While the figure axis draws a locus like this as the average for a long interval, the 
shape changes from year to year and, other than that, there exist many minor variations 
with shorter periods. Those minor variations are considered to be significant but not 
due to the errors in the observations or in the process of the analysis, and therefore 
they are necessary to be studied more in detail. 

Finally, a simple model with the seasonal flow of the matter in the Earth has 
been introduced to explain the variation of the figure axis obtained above. Applying 
the same model also to the variation of the rotational speed, it is confirmed that the 
variation estimated is accordant with the one which is so far known, thus suggesting 
that the variations of the figure axis and of the rotational speed are brought about by 
the same cause as far as the seasonal variation is concerned. 
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